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Abstract: The vegetation ecosystem of the Qinghai—Tibet Plateau in China, considered to be the "natural 
laboratory" of climate change in the world, has undergone profound changes under the stress of global 
change. Herein, we analyzed and discussed the spatial-temporal change patterns and the driving 
mechanisms of net primary productivity (NPP) in the Qinghai—Tibet Plateau from 2000 to 2015 based on 
the gravity center and correlation coefficient models. Subsequently, we quantitatively distinguished the 
relative effects of climate change (such as precipitation, temperature and evapotranspiration) and human 
activities (such as grazing and ecological construction) on the NPP changes using scenario analysis and 
Miami model based on the MOD17A3 and meteorological data. The average annual NPP in the Qinghai— 
Tibet Plateau showed a decreasing trend from the southeast to the northwest during 2000—2015. With 
respect to the inter-annual changes, the average annual NPP exhibited a fluctuating upward trend from 
2000 to 2015, with a steep increase observed in 2005 and a high fluctuation observed from 2005 to 2015. 
In the Qinghai—Tibet Plateau, the regions with the increase in NPP (change rate higher than 10%) were 
mainly concentrated in the Three-River Source Region, the northern Hengduan Mountains, the middle 
and lower reaches of the Yarlung Zangbo River, and the eastern parts of the North Tibet Plateau, 
whereas the regions with the decrease in NPP (change rate lower than —10%) were mainly concentrated in 
the upper reaches of the Yarlung Zangbo River and the Ali Plateau. The gravity center of NPP in the 
Qinghai—Tibet Plateau has moved southwestward during 2000-2015, indicating that the increment and 
growth rate of NPP in the southwestern part is greater than those of NPP in the northeastern part. 
Further, a significant correlation was observed between NPP and climate factors in the Qinghai—Tibet 
Plateau. The regions exhibiting a significant correlation between NPP and precipitation were mainly 
located in the central and eastern Qinghai—Tibet Plateau, and the regions exhibiting a significant 
correlation between NPP and temperature were mainly located in the southern and eastern Qinghai—Tibet 
Plateau. Furthermore, the relative effects of climate change and human activities on the NPP changes in 
the Qinghai—Tibet Plateau exhibited significant spatial differences in three types of zones, i.e., the climate 
change-dominant zone, the human activity-dominant zone, and the climate change and human activity 
interaction zone. These research results can provide theoretical and methodological supports to reveal the 
driving mechanisms of the regional ecosystems to the global change in the Qinghai—Tibet Plateau. 
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1 Introduction 


Net primary productivity (NPP) denotes the total amount of organic dry matter accumulated by 
green plants with respect to the unit time and area (Bhaskar, 2000; Rachhpal et al., 2007; Jiao et 
al., 2018). NPP is the residual part of the total amount of organic matter produced by 
photosynthesis after the elimination of autotrophic respiration. As the core content of the carbon 
budget and climate change research, NPP serves as an important indicator to determine the 
ecosystem health and sustainable development level (Alexander et al., 2007; Zhang et al., 2014; 
Leandro et al., 2019). Therefore, analyzing and discussing the spatial-temporal change patterns 
and the driving mechanisms of NPP in a regional vegetation ecosystem not only are of 
significance to understand the response of this ecosystem to climate change, but also can provide 
a scientific basis for ecological environment protection, plant yield estimation, resource 
exploitation, effective management of natural resources, and formulation of the corresponding 
social and economic development strategies (Sun et al., 2014). 

Currently, NPP research is conducted abroad using climate models (such as the 
BIOME-BioGeoChemical Cycle (BIOME-BGC) model, Carnegie-Ames-Stanford Approach 
(CASA) model and Global Production Efficiency (GLO-PEM) model) and remote sensing 
technology, whereas domestic scholars have utilized foreign empirical models (such as the Ball— 
Berry model and Boreal Ecosystem Productivity Simulator (BEPS) model) to comprehensively 
study the NPP in different spatial scale zones of China, such as northwest region, northeast region 
and southwest region (Mao et al., 2015; Guo et al., 2017; Han et al., 2018; Jiao et al., 2018; Xu et 
al., 2019; Zhu et al., 2019). Furthermore, the driving mechanisms on the response of NPP to 
climate change differ because of the differences in data sources, research areas and 
spatial-temporal scales (Chandrasekaran and Swamy, 2002; Travis et al., 2015). NPP globally 
increased from 1980 to 2000; however, the reasons varied in different regions (Li et al., 2018; Ji 
et al., 2019). Under the background of global change, climate in China has undergone significant 
changes, and the NPP of terrestrial vegetation exhibited an increasing trend (Chirici et al., 2007; 
Gao et al., 2018; Wang and Li, 2018; Wang et al., 2018). Climate change weakens the stress of 
hydrothermal conditions on vegetation growth. However, the response of vegetation NPP in 
China to climate change shows obvious regional differences (Guo et al., 2017; Li et al., 2017; 
Wang et al., 2017; Wang et al., 2018; Fatemeh et al., 2019). Therefore, further studies are required 
to investigate the driving mechanisms of NPP under the stress of global change in different 
regions (Edward et al., 2013). 

The unique geographical location, altitude, climate and geomorphological conditions of the 
Qinghai—Tibet Plateau in China result in various climatic and ecosystem types. As a special 
geographic unit, the Qinghai—Tibet Plateau exhibits an important role in regulating and indicating 
global change; it is considered to be the "natural laboratory" for studies of climate change in the 
world (Haberl et al., 2004; Huang et al., 2008; Sun et al., 2014; Wang et al., 2016). Although 
some previous studies have used conventional analysis methods (such as the correlation analysis 
model) to investigate the NPP change patterns in the Qinghai—Tibet Plateau, the relationships 
between NPP and climate factors and the relative effects of climate change and human activities 
on the NPP changes remain unclear (Sun et al., 2014; Huang et al., 2018; Jeremy et al., 2018; 
Sunil et al., 2019). Further, the "gravity center" concept, which indicates the point at which the 
force of gravity present on each part of an object exerts a resultant force, is derived from physics 
(Liu et al., 2011). However, in geography, the "gravity center" concept has been extensively 
applied to population, economy, food, land use, ecological environment, regional development 
and other fields (Liu et al., 2011). The spatial variation of the gravity center can better reflect the 
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variation degree and trend of geographical phenomena. 

Therefore, in this study, a gravity center model was introduced to analyze the spatial and 
temporal change patterns and the driving mechanisms of NPP in the Qinghai—Tibet Plateau from 
2000 to 2015 in combination with the conventional analysis methods. Furthermore, the relative 
effects of climate change and human activities on the NPP changes were also quantitatively 
distinguished. This study provides data and theoretical support to reveal the evolution patterns of 
the ecosystems in the Qinghai—Tibet Plateau and the response mechanisms of the ecosystems to 
global change. 


2 Study area and data sources 


2.1 Study area 


The Qinghai—Tibet Plateau extends from the Pamir Plateau and Karakorum Mountains in the west 
to the Hengduan Mountains in the east and from the Himalaya Mountains in the south to the 
Kunlun Mountains, Altun Mountains and Qilian Mountains in the north (26°00'12'—39°46'50’'N, 
73°18'52"-104°46'59”E). The plateau exhibits a total area of approximately 2.5x107 km? and an 
average altitude of 4000 m or more and is referred to as "the third pole". The annual mean 
temperature is lower than 0°C in the hinterland of the plateau; furthermore, the temperature 
decreases with increasing altitude and latitude. The mean annual precipitation decreases from 
2000 mm in the southeast to 50 mm in the hinterland of the plateau (Sun et al., 2014). In summer, 
abundant precipitation can be observed in most regions of the southeastern plateau because of the 
influence of the monsoon and humid airflow from the southwest Indian Ocean. In winter, strong 
winds (lasting for up to 200 d in some places) and scarce precipitation can be observed because of 
the influence of the dry westerlies in the plateau. Further, complex and diverse vegetation 
landscape types, including forests, meadows, grasslands and deserts, appear successively in the 
horizontal belt spectrum of the Qinghai—Tibet Plateau, and the climate gradually changes from the 
marine wet type in the southeast to the continental arid type in the hinterland of the plateau (Yang 
et al., 2015). 


2.2 Data sources 


The NPP during 2000-2015 obtained from the MODI17A3 dataset is available at 
https://Ipdaac.usgs.gov/. The spatial and temporal resolutions were 1 km and 1 a, respectively. 
The data (including the daily precipitation and mean daily temperature) of 229 meteorological 
stations in the Qinghai—Tibet Plateau and its surrounding areas, were obtained from the China 
Meteorological Data Network (http://data.cma.cn/). The MOD17A3 data were mosaicked and 
reprojected using the MODIS Projection Tool. Further, the annual mean temperature and annual 
precipitation were interpolated into grids with a spatial resolution of 1 kmx1 km using the 
Cokriging method in ArcGIS 10.3. Considering the uneven spatial distributions of the 
metrological stations in the western Qinghai—Tibet Plateau, we constructed the interpolated 
stations for the entire Qinghai—Tibet Plateau based on the China Meteorological Forcing Dataset 
(http://westdc.westgis.ac.cn), especially in the central and western regions. Locations of the 
meteorological and interpolated stations in the Qinghai—Tibet Plateau and its surrounding areas 
are shown in Figure 1. 


3 Methods 


3.1 Change trend analysis 


Linear regression analysis can be used to eliminate the influence of extreme climate in specific 
years to some extent (Tian et al., 2019). Therefore, this study adopted linear regression analysis to 
analyze the time series change trend of NPP, which can be calculated as follows: 
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Fig. 1 Spatial distributions of the average annual NDVI (during 2000-2015) in the Qinghai—Tibet Plateau and 
locations of the meteorological stations and interpolated stations in the Qinghai—Tibet Plateau and the 
surrounding areas 
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where slope denotes the gradient of change trend for vegetation NPP; n is the number of years or 
the whole study period (16 a); i is a variable of time and can be given as 7=1, 2, 3, ..., 16; and 
NPP; is the NPP in the i® year. The change rate of NPP can be given as follows: 
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where NPPer refers to the change rate of NPP in percentage (%); NPPmean is the average annual 
NPP during the period 2000-2015; and n is the number of years. 


3.2 Correlation analysis 


Correlation analysis is mainly used to indicate the correlation degree and direction between 
variables (Huang et al., 2008; Liu et al., 2011; Li and Zhao, 2013; Wang et al., 2016). In this 
study, we analyzed the correlation and significance level between NPP and climate factors 
(temperature and precipitation) using the Pearson's correlation coefficient based on pixels at 
confidence levels of 0.05 and 0.01. The correlation coefficient, partial correlation coefficient and 
complex correlation coefficient were calculated as Equations 3—5, respectively: 
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where rx, is the correlation coefficient of variables i and j; x;and y; denote the values of variables 


x2 


x and y in the i year, respectively; and x and y denote the average values of variables i and j, 


respectively. The value of r ranges from —1 to 1. A value greater than 0 indicates a positive 
correlation between variables i and y, while a value lower than 0 indicates a negative correlation. 
In Equations 4 and 5, rxy1,2 refers to the partial correlation coefficient of x and y after fixing y2; 
rxyı refers to the correlation coefficient of variables x and yı; rx;2 refers to the correlation 
coefficient of variables x and y2; ry1,2 refers to the correlation coefficient of variables yı and y2; 
Rxy1y2 refers to the complex correlation coefficient of variables x, yı and y2; and rxy2,1 refers to the 
partial correlation coefficient of variables x and y2 after fixing yı. In addition, the T test and F test 
were used to analyze the significance of the partial correlation coefficient and complex 
correlation coefficient. The formula can be given as follows: 


T =— 22 __xJn—m-1, (6) 
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where n denotes the number of samples and m denotes the number of variables. 
3.3 Gravity center model 


In a gravity center model, z; is considered as the attribute value of the number i plane space unit 
(grid) (Liu et al., 2011). When the Cartesian coordinate is set as (x; yi), the spatial mean of the 


region comprising three plane space units (grids) is defined as a Cartesian coordinate point (x, y) 
and can be expressed as: 


y=, (8) 


y=5 (9) 


With respect to this spatial phenomenon, an unbalanced distribution or "deviation of the gravity 
center" can be observed when the regional gravity center of a spatial phenomenon significantly 
differs from the geometric center. The deviation direction of the gravity center is toward the 
"high-density" area of the spatial phenomenon, and the deviation distance of the gravity center 
can be used to describe the degree of unbalanced distribution. Further, the mitigation trajectory of 
the gravity center in different periods can reflect the evolution process of the spatial phenomenon 
distribution. The formulas for obtaining the mitigation direction and mitigation distance of the 
gravity center are given as follows: 


0= arctan( %2 >t +m, (10) 


Xim ~ x, 


d,, = pS) + Oan IY 2 (11) 


where 0 and dm refer to the mitigation direction and mitigation distance of the gravity center, 
respectively; Ym and y, refer to the latitude coordinates of the gravity center in time t+m and time 
t, respectively; and x;+m and x; refer to the longitude coordinates of the gravity center in time +m 


and time ¢, respectively. 
3.4 Miami model 


In this study, the Miami model was adopted to estimate the potential vegetation NPP. This model 
based on Liebig's minimum method was the first NPP estimation model and has been widely 
adopted at home and abroad. It includes a function of the annual mean temperature and annual 
precipitation (Mao et al., 2015). 
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where NPPetimate is the potential vegetation NPP; T refers to the annual mean temperature (°C); 
and P refers to the annual precipitation (mm). 


3.5 Method for distinguishing the relative effects of climate change and human activities 
on the NPP changes 


Climate change and human activities are the main factors that affected the dynamic changes in 
vegetation NPP. According to the study of Mao et al. (2015), the relative effects of climate 
change and human activities on the NPP changes could be quantitatively evaluated and 
distinguished. In this study, three types of NPP changes were identified and obtained. The first 
type (actual NPP (NPPactual)), derived from MOD17A3, was based on the combined effects of 
climate change and human activities; the second type (potential NPP (NPPclimate)), obtained from 
the Miami model, was based on the effects of climate change; and the third type (human 
activity-related NPP (NPPhuman)) was the residual of the NPPelimate and NPPactual, which can be 
mainly attributed to human activities such as grazing and ecological construction. 

Further, we calculated the change trends of NPPeiimate and NPPhuman (Ke and Kn, respectively) to 
assess and distinguish the relative effects of climate change and human activities on the NPP 
changes. When K. is positive, climate change is conducive to vegetation restoration. In contrast, 
climate change exacerbates vegetation degradation when K. is negative. Human activities 
exacerbate vegetation degradation when Kp is positive, whereas human activities are conducive to 
vegetation restoration when Kp is negative. The change trend of NPPactual (Ka) could indicate 
either vegetation restoration or degradation. When Ka is positive, the vegetation is observed to 
undergo restoration, whereas the vegetation is observed to undergo a degradation process when 
Ka is negative. As shown in Table 1, six possible scenarios for the NPP changes have been 
defined (Sun et al., 2014; Wang et al., 2018). 


Table 1 Relative effects of climate change and human activities on the NPP changes 


p Relative effects of Relative effects of human 
Change trend of NPP Scenario Ke Ka climate change (%) activities (%) 
Scenario 1 >0 >0 100 0 
K0 Scenario 2 <0 <0 0 100 
(vegetation restoration) 
Scenario 3 >0 <0 Combined effects Combined effects 
Scenario 4 >0 >0 0 100 
K,<0 : 
(vegetation degradation) Scenarios = E 109 9 
Scenario 6 <0 >0 Combined effects Combined effects 


Note: K, refers to the change trend of actual NPP (NPPactuai); and K. and Ka refer to the change trends of potential NPP (NPPetimate) and 
human activity-related NPP (NPPhuman), respectively. 


4 Results 


4.1 Spatial variation characteristics of NPP 


4.1.1 Spatial distribution patterns of average annual NPP during 2000-2015 

Figure 2 denotes that the spatial distribution of average annual NPP significantly differed in the 
Qinghai—Tibet Plateau during 2000-2015. Generally, a decreasing trend can be observed from the 
southeast to the northwest; this was related to the spatial and temporal distributions of the 
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regional hydrothermal and climate conditions. During the study period, the average annual NPP in 
the Qinghai—Tibet Plateau was 167.52 g C/(m?-a). Further, the average annual NPP for majority of 
the study area, such as the North Tibet Plateau, the middle and upper reaches of the Yarlung 
Zangbo River, and the Ali Plateau, was lower than 400.00 g C/(m*a). A high-value zone of 
average annual NPP can be observed in the southeastern part of the Qinghai—Tibet Plateau. The 
regions exhibiting an average annual NPP of greater than 600.00 g C/(m?.a) were mainly 
distributed in the southern part of the Hengduan Mountains. The regions with an average annual 
NPP of 200.00-400.00 g C/(m?-a) were mainly distributed in the West Sichuan Plateau. A 
decreasing trend for the average annual NPP was observed from the south to the north and from 
the west to the east in the Hengduan Mountains. The regions exhibiting high average annual NPP 
values in the Hengduan Mountains were mainly concentrated in the southern part, which can be 
attributed to the vegetation type, topography and climate in this region. 
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Fig. 2 Spatial distributions of average annual NPP (net primary productivity) in the Qinghai—Tibet Plateau 
during 2000-2015 


4.1.2 Spatial change patterns of average annual NPP 


The change rate of NPP in the Qinghai—Tibet Plateau from 2000 to 2015 was analyzed. During 
the study period, the regions exhibiting the change rate of NPP in the range of 10%-—30% in the 
Qinghai—Tibet Plateau were the largest, accounting for 45.55% of the whole study area, followed 
by the regions exhibiting the change rate of NPP in the range from —10% to 10%, accounting for 
33.39% (Fig. 3). The regions with the change rate of NPP lower than —30% were the smallest, 
only accounting for 0.91%. These results denote an increasing trend with respect to the 
inter-annual changes of NPP in the Qinghai—Tibet Plateau from 2000 to 2015, indicated by the 
stable (-10%<change rate of NPP<10%) and increasing (change rate of NPP>10%) regions, 
accounting for 94.20% of the total study area. 

Figure 4 shows that the stable regions with the change rate of NPP in the range from —10% to 
10% were mainly distributed in the eastern and southern parts of the Qinghai—Tibet Plateau, 
including the northern Hengduan Mountains and the southern Ali Plateau. The regions exhibiting 
the change rate of NPP greater than 30% were mainly distributed in the eastern part of the North 
Tibet Plateau, the lower reaches of the Yarlung Zangbo River and the middle of the Hengduan 
Mountains. The regions exhibiting the change rate of NPP in the range of 10%—30% were the 
most widely distributed and mainly concentrated in the Three-River Source Region, the northern 
Hengduan Mountains, the middle and lower reaches of the Yarlung Zangbo River, and the middle 
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and eastern parts of the North Tibet Plateau. The regions with NPP reduction (change rate of NPP 
lower than —10%) were mainly distributed in the upper reaches of the Yarlung Zangbo River and 
the Ali Plateau. 
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Fig. 3 Area percentages of regions with different grades in the change rate of NPP in the Qinghai—Tibet Plateau 
during 2000-2015 
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Fig. 4 Spatial distributions of regions with different grades in the change rate of NPP in the Qinghai-Tibet 
Plateau during 2000-2015 


The gravity center could effectively denote the imbalance and bias of NPP changes with respect 
to the spatial distribution. In this study, the gravity centers of NPP for each year and average NPP 
of the total study period (16 a) were calculated (Fig. 5). It can be seen that the gravity centers of 
NPP were mainly concentrated in the southeastern part of the Qinghai-Tibet Plateau. This 
indicated that the total NPP in the southeastern part of the Qinghai—Tibet Plateau was larger when 
compared with that in the northwestern part. This result was consistent with the conclusion 
obtained from Figure 2 that the average annual NPP decreased from the southeast to the northwest 
in the Qinghai—Tibet Plateau. 

In this study, we established the polar coordinate system by adopting the mean NPP gravity 
center (during 2000-2015) in the study area as the origin (representing the average state) and 
using the migration distance and migration angle of the NPP gravity center to the origin as the 
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polar radius and polar angle, respectively, to further analyze the spatial distribution patterns of the 
gravity centers for NPP in the Qinghai—Tibet Plateau. The number of the NPP gravity centers in 
different quadrants of the coordinate system and its percentage were estimated. The results shown 
in Figure 6 denoted that 37.50% of the NPP gravity centers were located in the southwest 
quadrant, which were followed by those in the northeast quadrant with a proportion of 25.00%. 
This indicated that the increment and growth rate of NPP in the southwestern and northeastern 
parts of the Qinghai—Tibet Plateau were relatively large during the study period. In addition, the 
sum of the gravity center number in the two aforementioned quadrants accounted for 62.50% of 
the total gravity center number. Further, the proportions of the NPP gravity centers in the 
northwest and southeast quadrants were relatively small at 18.75%. Thus, we divided the study 
area into two parts based on the coordinate axis. The gravity center number of NPP in the 
northern half quadrant (43.75%) was slightly smaller than that in the southern half quadrant 
(56.25%). This indicated that the increment and growth rate of NPP in the southern part of the 
Qinghai—Tibet Plateau were larger than those in the northern part during 2000-2015. Meanwhile, 
the gravity centers of NPP had larger migration distances when compared with the origin in 2002, 
2003, 2005 and 2014, which indicated a significant increment of NPP in the northeastern and 
southwestern parts of the Qinghai—Tibet Plateau. In 2006, 2007 and 2008, the mitigation distances 
of the gravity center of NPP were smaller, implying that the changes in NPP in each parts of the 
Qinghai—Tibet Plateau were relatively consistent in these three years. 
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Fig. 5 Distributions of the gravity centers of NPP in the Qinghai—Tibet Plateau during 2000-2015 


Further, we analyzed the migration directions of the NPP gravity center at different temporal 
scales during 2000-2015 to explore the spatial variations of the gravity center of NPP in the 
Qinghai—Tibet Plateau. Considering the influences of extreme climate and other factors on the 
annual change in vegetation NPP, we adopted three temporal scales (3, 5 and 16 a) to analyze the 
migration directions of the gravity center of NPP in the Qinghai—Tibet Plateau. Figure 7 denotes 
that the gravity center of NPP moved eastward from 2000-2002 to 2003-2005, indicating that the 
increment and growth rate of NPP in the eastern part of the study area were higher than those in 
the western part. The gravity center of NPP during 2006-2008 moved toward the southwest 
direction when compared with that during 2003-2005, indicating that the growth rate of NPP in 
the southwestern part was greater than that in the northeastern part. Further, when compared with 
the gravity center of NPP during 2006-2008, the gravity centers during 2009-2011 and 2012- 
2015 moved westward, suggesting that the increment of NPP in the western part of the Qinghai- 
Tibet Plateau was larger than that in the eastern part during 2009-2015. On a 5-a temporal scale, 
the gravity center of NPP in the Qinghai—Tibet Plateau initially moved to the southeast direction 
from 2000-2004 to 2005-2009 and subsequently moved to the west direction from 2005-2009 to 
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2010-2015. Thus, in the Qinghai—Tibet Plateau, the increments of NPP in the southeastern part 
during 2005-2009 and in the western part during 2010-2015 were greater than those in the 
northwestern part during 2000-2004 and in the eastern part during 2005-2009, respectively. 
Generally, during the period 2000-2015, the gravity center of NPP in the Qinghai—Tibet Plateau 
moved southwestward, indicating that the increment and growth rate of NPP in the southwestern 
part were larger than those in the northeastern part. 
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Fig. 6 Polar coordinates of the gravity centers of NPP in the Qinghai—Tibet Plateau during 2000-2015. The 
green dots indicate the gravity centers. 
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Fig. 7 Migration trajectory of the gravity centers of NPP in the Qinghai—Tibet Plateau during 2000-2015 at 
different temporal scales (3-, 5- and 16-a temporal scales) 


4.2 Inter-annual variations of average annual NPP 


From 2000 to 2015, the average annual NPP in the Qinghai—Tibet Plateau showed an increasing 
trend with a slight fluctuation (Fig. 8). The fluctuation range was 0.00-32.88 g C/(m?.a). During 
the study period, the maximum NPP value was 182.95 g C/(m2.a) in 2013, followed by 180.88 g 
C/(m?.a) in 2006, whereas the minimum value was 150.07 g C/(m*a) in 2000. Before 2005, the 
average annual NPP in the Qinghai—Tibet Plateau was less than 160.00 g C/(m?.a) and the 
increase trend was minimal. After 2005, however, the average annual NPP increased significantly, 
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with the mean value of 170.00 g C/(m?.a). In 2006 and 2013, the average annual NPP values were 
180.88 and 182.95 g C/(m.a), respectively. This increasing trend of average annual NPP mainly 
resulted from the implementation of the ecological projects of "returning farmland to grassland" 
and "returning grazing land to grassland". 
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Fig. 8 Inter-annual variations of average annual NPP in the Qinghai—Tibet Plateau during 2000-2015 


4.3 Correlation analyses of NPP with precipitation and temperature 
4.3.1 Changes of climate factors during 2000-2015 


As shown in Figure 9, the annual precipitation of the Qinghai—Tibet Plateau during 2000-2015 
ranged from 379.36 to 541.10 mm. In addition, the average annual precipitation was 490.65 mm 
during the study period, with the maximum value observed in 2010 and the minimum value in 
2007. During the study period, a slight decrease was observed in the annual precipitation in the 
Qinghai—Tibet Plateau; further, the average annual precipitation decreased from the southeast to 
the northwest with respect to the spatial distribution. During 2000-2015, the annual mean 
temperature significantly increased in the Qinghai—Tibet Plateau, ranging from 0.67°C (2002) to 
2.98°C (2007). Meanwhile, the mean annual temperature during the study period was 1.48°C. 
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Fig. 9 Inter-annual variations of (a) annual precipitation and (b) annual mean temperature in the Qinghai—Tibet 
Plateau during 2000-2015 


4.3.2 Correlation analyses of NPP with precipitation and temperature 


Because the influences of precipitation and temperature on the NPP are often intersectional and 
dependent, the partial correlation coefficient and complex correlation coefficient between NPP 
and precipitation and between NPP and temperature in the Qinghai—Tibet Plateau during 2000- 
2015 were calculated pixel by pixel (Fig. 10). 

As shown in Figure 10a, the partial correlation coefficient between NPP and annual 
precipitation ranged from —0.86 to 0.91. The regions exhibiting positive partial correlation 
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coefficients between NPP and annual precipitation accounted for 54.57% of the total study area 
and were mainly distributed in the central and northeastern of the Qinghai—Tibet Plateau. In 
contrast, the regions with negative partial correlation coefficients between NPP and annual 
precipitation were mainly concentrated in the southern part of the North Tibet Plateau, the upper 
reaches of the Yarlung Zangbo River Basin and the Ali Plateau. 
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Fig. 10 Partial correlation coefficient, multiple correlation coefficient and significant test results of NPP with 
precipitation and temperature in the Qinghai—Tibet Plateau. (a), partial correlation coefficients between NPP and 
precipitation; (b), significance of T test for the partial correlation coefficients between NPP and precipitation; (c), 
partial correlation coefficients between NPP and temperature; (d), significance of T test for the partial correlation 
coefficients between NPP and temperature; (e), complex correlation coefficients between NPP and climate factors 
(precipitation and temperature); (f), significance of F test for the complex correlation coefficients between NPP 
and climate factors. Coeff., correlation coefficient. 


As shown in Figure 10c, the regions with positive partial correlation coefficients between NPP 
and temperature accounted for 47.13% of the total study area and were mainly distributed in the 
eastern and southern parts of the Qinghai—Tibet Plateau. Meanwhile, the regions with negative 
partial correlation coefficients accounted for 52.87% of the total study area and were mainly 
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distributed in the western part of the Three-River Source Region, the upper reaches of the Yarlung 
Zangbo River Basin, the southern part of the Hengduan Mountains and its surrounding areas. In 
addition, the significances of partial correlation coefficients between NPP and precipitation and 
between NPP and temperature were evaluated using the T tests. 

The results shown in Figures. 10b and d revealed that for the partial correlation coefficients 
between NPP and precipitation, 15.83% of the total study area passed the T test at P<0.05 level, 
whereas 5.34% of the whole area passed the T test at P<0.01 level. Thus, the regions exhibiting 
the extremely significant correlation between NPP and precipitation were mainly distributed in 
the western and eastern parts of the Three-River Source Region. For the partial correlation 
coefficients between NPP and temperature, approximately 12.65% of the total study area passed 
the T test at P<0.05 level, whereas 6.13% of the total area passed the T test at P<0.01 level. The 
regions with the extremely significant correlation between NPP and temperature were mainly 
distributed in the lower reaches of the Yarlung Zangbo River Basin. 

As depicted in Figure 10e, with respect to the complex correlation coefficients between NPP and 
climate factors (precipitation and temperature), the coefficient values ranged from 0.00 to 0.94. In 
the significance of F test for the complex correlation coefficients between NPP and climate factors 
(Fig. 10f), approximately 16.32% of the total area passed the F test at P<0.05 level, whereas 
5.89% of the whole study area passed the F test at P<0.01 level. Therefore, the regions exhibiting 
extremely significant complex correlations between NPP and climate factors were mainly 
distributed in the middle and western parts of the Three-River Source Region, the upper and lower 
reaches of the Yarlung Zangbo River Basin, and the West Sichuan Plateau. 


4.4 Effects of climate change and human activities on the NPP changes 


As shown in Figure 11, the regions with vegetation restoration (increase of NPP) accounted for 
84.54% of the Qinghai—Tibet Plateau during 2000-2015, whereas the regions with vegetation 
degradation (decrease of NPP) occupied 15.46% of the whole study area. The regions with 
vegetation restoration caused by climate change accounted for 46.01% of the total study area and 
were mainly distributed in the western and central parts of the North Tibet Plateau and the 
Three-River Source Region. The regions with vegetation restoration caused by human activities 
accounted for 34.33% of the total study area and were mainly concentrated in the middle and 
lower reaches of the Yarlung Zangbo River Basin and the Hengduan Mountains. Furthermore, the 
combined effects of climate change and human activities resulted in a small area of vegetation 
restoration (3.81% of the study area), which was mainly distributed in the southern part of the 
Three-River Source Region. Meanwhile, the regions with vegetation degradation caused by 
climate change accounted for 5.05% of the total study area and were discontinuously distributed 
in the upper reaches of the Yarlung Zangbo River Basin and the northern part of the Hengduan 
Mountains. The regions with vegetation degradation caused by human activities accounted for 
10.20% of the study area and were mainly concentrated in the Ali Plateau. Moreover, the regions 
with vegetation degradation that can be attributed to climate change and human activities were 
mainly distributed in the northern margin of the Yarlung Zangbo River Basin. 

Thus, as shown in Figure 11, climate change was considered to be the dominant influencing 
factor in the Al and A2 regions (climate change-dominant zone), the interaction zone (B) was 
considerably influenced by the combined effects of climate change and human activities, and the 
C1 and C2 regions (human activity-dominant zone) were dominated by human activities. 


5 Discussion 


5.1 Spatial and temporal distribution patterns of NPP 


The spatial distribution of NPP in the Qinghai—Tibet Plateau showed a decreasing trend from the 
southeast to the northwest during 2000-2015, which was mainly influenced by the spatial 
distribution patterns of precipitation and temperature. Because of the Pacific Monsoon, the 
Hengduan Mountains in the southeastern part of the Qinghai—Tibet Plateau exhibit abundant 
precipitation, diverse vegetation types, broad-leaved forests, coniferous forests and shrubs, and 
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high vegetation NPP. The lower reaches of the Yarlung Zangbo River are affected by the warm 
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Fig. 11 Distinction of the relative effects of climate change and human activities on the NPP changes (increase 
or decrease). Al and A2 refer to the regions mainly affected by climate change; B refers to the region affected by 
the combined effects of climate change and human activities; C1 and C2 refer to the regions mainly affected by 
human activities; L1, L2, L3 and L4 refer to the boundary lines for different regions. 


and wet airflow in the Indian Ocean. The forest cover is high and the NPP is high in rainforests, 
monsoon forests and broad-leaved forests. The NPP in most areas of the northern Tibet Plateau is 
relatively low because of the high altitude, low temperature, low precipitation and numerous 
vegetation types, such as grasslands, alpine meadows and sparse shrubs (Sun et al., 2014). Deserts 
and desert grasslands with the lowest NPP are widespread in the Ali Plateau because of the 
scarcity of annual precipitation and low temperature. During 2000-2004, the average annual NPP 
in the Qinghai—Tibet Plateau has shown an upward trend. After 2004, the average annual NPP has 
shown a steep increase. This increase can be attributed to the fact that the forest lands and 
grasslands have been considerably restored after the implementation of the projects of "returning 
farmland to grassland" and "returning grazing land to grassland". The average annual NPP in the 
Qinghai—Tibet Plateau presented a sharp increase in 2005, with high fluctuations observed during 
2005-2015. 


5.2 Effects of climate factors on NPP 


Based on the comprehensive analyses of NPP, precipitation and temperature, we can infer that the 
correlation between NPP and climate factors is significant in different regions and may be related 
to the altitude, climate type and vegetation type. In the central, northeastern and western parts of 
the Qinghai—Tibet Plateau, NPP was positively correlated with precipitation because the altitude 
of these regions was higher than 4000 m a.s.l. and the precipitation was low. Therefore, 
precipitation became the dominant factor that restricts the growth of vegetation. In addition, the 
grasslands and alpine meadows in these regions demonstrated a high sensitivity to precipitation 
(Yang et al., 2015). However, NPP was negatively correlated with precipitation in the 
southeastern part of the Hengduan Mountains, and the middle and lower reaches of the Yarlung 
Zangbo River Basin. Precipitation was abundant in these regions, so it was not the main factor 
that restricts the growth of vegetation. Because of the high levels of precipitation in these regions, 
the rainfall erosivity was high, the slope was steep and the soil erosion phenomena (freeze-thaw 
erosion and hydraulic erosion) were serious. In addition, the negative correlations between NPP 
and precipitation in some regions (especially in the Hengduan Mountains) were associated with 
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lower solar radiation due to increased precipitation (Liu et al., 2011). Thus, the regional 
vegetation ecosystem will be destroyed to a certain extent. 

The regions with positive correlations between NPP and temperature were mainly distributed in 
the southern part of the Tibet Plateau and the northern part of the Hengduan Mountains. These 
regions exhibited high altitude, abundant precipitation and low temperature (restraining the 
growth of vegetation); therefore, the higher the temperature, the better the growth of vegetation 
(strong photosynthesis). The regions that exhibited negative correlations between NPP and 
temperature were mainly concentrated in the central and northern parts of the Qinghai-Tibet 
Plateau because of scarce precipitation in these regions. Although the increase in temperature was 
conducive to the plant photosynthesis, regional drought was exacerbated to a certain extent, 
inhibiting vegetation growth. 


5.3 Relative effects of climate change and human activities on the NPP changes 


Climate change and human activities were the main driving forces of the NPP changes in the 
Qinghai-Tibet Plateau. In this study, the NPPactuaı (MOD17A3) and NPPetimate (reversed by Miami 
model) were used to establish public indicators to unify the effects of climate change and human 
activities on the vegetation degradation or restoration to a comparable level. The results showed 
that the relative effects of climate change and human activities on the NPP changes in the 
Qinghai—Tibet Plateau exhibited significant spatial differences. The study area can be divided into 
three types of zones, i.e., the climate change-dominant zone, the human activity-dominant zone, 
and the climate change and human activity interaction zone. The climate change-dominant zone 
could be divided into two parts: Al (a region with the increase in NPP induced by climate change) 
and A2 (a region with the decrease in NPP induced by climate change). Figure 12 shows the 
spatial distributions of change trend coefficients of precipitation and temperature in the Qinghai— 
Tibet Plateau during 2000-2015. 
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Fig. 12 Spatial distributions of change trend coefficients of (a) precipitation and (b) temperature in the Qinghai— 
Tibet Plateau during 2000-2015 


In the Al region, precipitation and temperature exhibited a significant increasing trend (Fig. 
12) and water resources were abundant, contributing to the photosynthesis of vegetation and 
considerably promoting the growth of vegetation (Huang et al., 2008). In the A2 region, however, 
a significant increase in temperature was observed due to the scarcity of precipitation, resulting in 
different degrees of drought and thus destroying the regional vegetation ecosystem. The human 
activity-dominant zone can also be divided into two parts: C1 (a region with the decrease in NPP 
induced by human activities) and C2 (a region with the increase in NPP induced by human 
activities). Grazing and farming were the main human activities in these regions. However, with 
the implementation of the ecological projects of "returning farmland to grassland" and "returning 
grazing land to grassland", the effects of human activities on the NPP changes in different regions 
showed different directions (inhibition or promotion). Overgrazing and reclamation played a 
leading role in the C1 region, whereas ecological protection policies played an essential role in 


the C2 region (Sun et al., 2014). In addition, in the Qinghai—Tibet Plateau, there existed an 
interaction zone (B) where the NPP changes were dominated by the combined effects of climate 
change and human activities. This region was a transitional zone in which human effects were 
gradually weakening while climate effects were gradually increasing. Thus, significant temporal 
and spatial differences can be observed in the relative effects of climate change and human 
activities on the NPP changes in the Qinghai—Tibet Plateau. 


6 Conclusions 


In this study, the spatial-temporal change patterns and the driving mechanism of NPP in the 
Qinghai—Tibet Plateau during 2000-2015 were discussed. From 2000 to 2015, the average annual 
NPP exhibited a decreasing trend from the southeast to the northwest over the whole Qinghai- 
Tibet Plateau. In terms of the inter-annual changes of NPP, a fluctuating upward trend was 
observed during the study period, with a steep rise observed in 2005 and a high fluctuation 
observed during 2005-2015. During 2000-2015, the increment and growth rate of NPP in the 
southwestern part were greater than those in the northeastern part. The regions that exhibited 
significant correlations between NPP and precipitation were mainly located in the central and 
eastern parts of the Qinghai—-Tibet Plateau, whereas the regions that exhibited significant 
correlations between NPP and temperature were mainly distributed in the southern and eastern 
parts. The relative effects of climate change and human activities on the NPP changes in the 
Qinghai—Tibet Plateau exhibited significant spatial differences. The study area can be divided into 
three types of zones, i.e., the climate change-dominant zone, the human activity-dominant zone, 
and the climate change and human activity interaction zone. 
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